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Abstract 

We study the Plancherel-Rotach asymptotics of four families of orthogonal poly- 
nomials, the Chen-Ismail polynomials, the Berg-Valent polynomials, the Conrad- 
Flajolet polynomials I and II. All these polynomials arise in indeterminate moment 
problems and three of them are birth and death process polynomials with cubic 
or quartic rates. We employ a difference equation asymptotic technique due to 
Wang and Wong. Our analysis leads to a conjecture about large degree behav- 
ior of polynomials orthogonal with respect to solutions of indeterminate moment 
problems. 
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1 Introduction 



The Plancherel-Rotach asymptotics for orthogonal polynomials refer to asymptotics of 
the orthogonal polynomials scaled by the largest zero. There are three types of Plancherel- 
Rotach asymptotics depending on whether we are in the oscillatory region between the 
largest and smallest zeros, in the exponential region beyond the largest (or smallest) 
zero, or in the transition region near the largest (or smallest) zero. If the extreme zero 
moves as the polynomial degree increases, the last type of asymptotics is related to the 
soft edge asymptotics in the random matrix theory and describes the eigenvalues of large 
random Hermitian matrices near the tail end, see [2]. The Plancherel-Rotach asymp- 
totics for Hermite and Laguerre polynomials are given in [20] and [2] . The corresponding 
asymptotics for the exponential weights were developed by using the Deift-Zhou steepest 
descent analysis for corresponding Riemann-Hilbert problems in [TU] and [H]. 

Every birth and death process leads to a family of orthogonal polynomials. The birth 
rates {A^} and the death rates {fin} are assumed to satisfy the conditions that A„ and 
fin+i are positive and /iq > 0. The recurrence relation of the birth and death process 
polynomials {Qn{x)} is 

- xQn{x) = \nQn+l{x) + HnQn-l{x) - i\n + IJ'n)Qn{x), n>0 (1.1) 

with the initial condition 

Qo{x) = 1, Qi(x) = (Ao + fio- x)/Xo. (1.2) 

It is known in the literature that some classical continuous and discrete orthogonal poly- 
nomials can be written as birth and death process polynomials with certain special rates. 
For example, one can find the birth and death rates for Laguerre polynomials, Jacobi 
polynomials and Charlier polynomials in [T3| Sec. 5.2]. Of course, for these classical 
polynomials, their weight functions are well-known and unique. However, for quite a few 
choices of A„ and the corresponding moment problem is indeterminate. This means 
that the orthogonality measures of such polynomials are not unique. For example. Berg 
and Valent [S] introduced polynomials associated with quartic rates 

A„ = (4n + l)(4n + 2)2(4n + 3), /i„ = (4n - l)(4n)2(4n + 1), (1.3) 

Conrad and Flajolet [SI E] considered the following two families with cubic rates 

Xn = {3n + c+l){3n + c + 2f, = (3n + c)2(3n + c + 1), c> 0, (1.4) 

A„ = (3n + c+ l)2(3n + c + 2), /i„ = (3n + c - l)(3n + c)^ c> 0; (1.5) 

see also [11] . The moment problems associated with the above rates are all indeterminate. 
In fact, the orthogonality measures are still unknown in the literature. 
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Chen and Ismail [7] considered the orthogonal polynomials rising from the following 
recurrence relation 

Tn+i{x) = xTn{x) - An^An" - 1) (1.6) 

with J-o(x) = 1 and J^i{x) = x. They indicated that these polynomials are also related 
to the indeterminate moment problems. Moreover, they found the following family of 
weight functions for J-'„(x) 

^ ^ v^r^[cos(pv/^) + cosh(pv/^)] ^ ^ ^ 

2[{cos(pya72) + cosh(pv/^)}2 - a2sin2(p^/;^) sinh2(pv/^)] ' 

(1.7) 

where a G (—1, 1) and p is a constant defined in (12.371) below; cf. [3 eq. (5.11)]. 

As the orthogonality measures are usually unknown for orthogonal polynomials re- 
lated to the indeterminate moment problems, it is difficult to apply the powerful Deift- 
Zhou steepest descent analysis for Riemann-Hilbert problems. Even for the Chen-Ismail 
polynomials, because the weight function (II. 7p is not analytic in the neighbourhood of 
M, it is difficult to apply the Riemann-Hilbert technique. So we intend to use the dif- 
ference equation method developed by Wong and his colleagues in [221 [211 [25] to get 
the asymptotic expansions. This will be the first time to apply this technique to derive 
Plancherel-Rotach asymptotics for orthogonal polynomials corresponding to indetermi- 
nate moment problems. It must be emphasized that we use the recurrence relation (or 
difference equations) only to achieve our asymptotic expansion in this paper. This is a 
significant improvement comparing with the previous works, where some extra asymp- 
totic information is needed to determine the asymptotic expansion; for example, see 
[I8l[2a[25]. 

Before we proceed to the next section, we would like to mention that one useful tool 
in estimating the largest and smallest zeros of orthogonal polynomials is the following 
theorem of Ismail and Li [T3] . 

Theorem 1.1. Let {Vn{x)} he a sequence of monic polynomials satisfying 

xVn{x) = Vn+lix) + anVn{x) + PnVn-lix), (1.8) 

with Pn > 0, for 1 < n < N and let {c„} be a chain sequence. Set 

B := max{a;„ : < n < A^}, and A := min{?/„ : < n < A^}, (1.9) 
where Xn and Vn, Xn > Vn, o-"^^ the roots of the equation 

{X - an){x - an-l)Cn = Pn, (1-10) 
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that is 



(a„ + an-i 



an - a, 



(1.11) 



Then the zeros ofVN{x) lie in {A,B). 

In practice, one usually choose the simple chain sequence c„ = 1/4 for all n in the 
above theorem. 

The present paper is organized as follows. In Section 2, we study the Chen-Ismail 
polynomials. After applying results in Wang and Wong [25] to obtain two linearly in- 
dependent solutions for the difference equation (11. 6p . we develop ideas used in Wang 
and Wong [23] to determine the coefficients of these two solutions. Then the Plancherel- 
Rotach asymptotic expansion is obtained, which is given in terms of Airy functions. 
Moreover, we study the limiting zero distribution as well as the behavior of the largest 
zero. In the next three sections, we follow the similar procedures as in Section 2 and 
derive asymptotic expansions for Berg-Valent polynomials and Conrad-Flajolet polyno- 
mials, respectively. Since Bessel-type asymptotic expansions appear in these cases, we 
need to apply the recent results by Cao and Li [6] to get the two linearly independent 
solutions of our difference equations. As Sections 2-5 are independent of each other, we 
shall use the same notation for the functions and variables but it will mean different 
things in different sections. We hope that this will not cause any confusion. In the last 
section of this paper, we list some remarks about the moment problem and formulate a 
conjecture about large degree (Plancherel-Rotach) behavior of orthogonal polynomials 
associated with indeterminate moment problems. In this section, we also apply asymp- 
totic techniques to show that the moment problem associated with the Conrad-Flajolet 
polynomials is indeterminate. 

2 Chen— Ismail polynomials 

Before we derive the Plancherel-Rotach asymptotics for J-'„(x), one can easily get a bound 
for the largest and smallest zeros from Theorem II. 1[ 

Proposition 2.1. Let Xn,k be zeros of J^n{x) such that Xn^i > 3^n,2 > ■ ■ ■ > a^n.n? then we 
have the following bounds for all n > 1 



Proof. Recall the recurrence relation of J^n(x) in (II. 6p and choose = 1/4 in Theorem 




and 




(2.1) 



11.11 Then the result follows. 



□ 
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2.1 Difference equation method 

To derive the Plancherel-Rotach asymptotics of J^n{x), we start from the recurrence 
relation in fll.6p and apply Wang and Wong's difference equation method developed in 
[211 [25]. To use their results, we transform the recurrence relation (11.61) into the following 
standard form they need 

Pn+l{x) - {AnX + Bn) Pn{x) + Pn-l{x) = 0. (2.2) 

Let 

K„ 2'^ r- (^1) r ('^) r (2.3) 



. 2 ^ ^ 
and J^n{x) = KnPn{x), then (II. 6p becomes 

Kn+lPn+l{x) = xKnPn{x) - An^{An'^ - l)Kn-lPn~l{x). (2.4) 

Since Kn+i = 4?7,^(4?7,^ — l)A'„„i, the above recurrence equation reduces to the standard 
form ([22D with v4„ = and 5„ = 0. As n — )■ oo, the recurrence coefficient satisfies 
the following expansion 

oo 

a 



n 

=0 



with 6 = 2 and 



11 . ^ 

-, a, = --. (2.6) 



Since An is of order 0(r;,~^) when n is large, to balance the term AnX in (12. 2p . we introduce 
X = u'^t and u = n + tq, where Tq is a constant to be determined. Then the characteristic 
equation for (12.21) is 

- aotX + 1 = (2.7) 
with ao given in (12. 6p . The roots of this equation are 



(2. 



and they coincide when the quantity inside the above square root vanishes, that is, 
t = t± = ±8. These points t± are called transition points for difference equations by 
Wang and Wong in [M[ [25] because in their neighbourhood the behaviors of solutions to 
(12. 2p changes dramatically. 

Since the polynomials P2n{x) and P2n+i{x) are even and odd functions, respectively, 
let us consider the asymptotics near the large transition point t+ = 8 only. According to 
the results in [23], we have the following proposition. 
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Proposition 2.2. When n is large, Pn{x) in fl2.2p can he expressed as 

Pn{x) = Ci(x)P„(x) + C2(x)(5„(x), 



(2.9) 



where Ci{x) andC2{x) are two n-independent functions. In the above formula, Pn{x) and 
Qn{x) are two linearly independent solutions of (12.21) satisfying the following Airy-type 
asymptotic expansions in the neighbourhood oft+ = 8 



Pn(^'t) 



and 



uu{t) 

t2 - 64 



uu{t)\ 
t2 - 64 y 



s=0 



V 6 



s=0 



(2,10) 



Bi(.fc/w)f:M^ + Bi'(.it/w)f;^"<''> 



U 6 



s=0 



U ' 6 



f2.ir 



where z/ = n + | anc? ?7(t) defined as 



,7r 



F(-,-l)-F(arcsinW-,-l^ 



log 



t + - 64 



-[-f/(t)]5=cos-^- 



V2t^ ,1 1 



B 



1-— ' 



2' 4' 



< t < 



Here F{(f, k^) is the elliptic integral of the first kind 

F{^,k^) = 

Bx{a,b) is the incomplete Beta function 



t > 8,(2.12) 
(2.13) 

(2.14) 



B,{a,b)= / y''-\l-yf~'dy 
Jo 

and the leading coefficients are given by 

Ao{U) = 1, MU) = 0. 



(2.15) 



(2.16) 



Proof. Recall the asymptotic expansion for An in (12. 5p . the equation (12. 2p falls into the 
case 6^0 and t^ ^ considered in [2l]. Following their approach, we choose 



ait+ 1 

Tn = = - and u = n + Tn. 

^ 26 2 



Then our proposition follows from the main Theorem in 



(2.17) 
□ 
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Remark 2.1. Note that the terms involving the eUiptic integral F{if, k"^) and the incom- 
plete Beta function Bx{a, b) in f l2.12p and f l2.13p can be written as 



F(|,-1)-F(arcsin -1) 



- 64 



,2 -,- = / ' ' ds. 



Moreover, one can verify that ?7(t) defined in fl2.12p and fl2.13p is a monotonically in- 
creasing function in the neighbourhood of 8. In fact, we have the following asymptotic 
formula 

[/(t) = + 0(t - 8)^ ast^8. (2.18) 
8 V 2 

2.2 Determination of Ci{x) and C2{x) 

In this subsection, we will determine the coefficients Ci{x) and C2{x) in Proposition 12.21 
via a matching method. To this end, we shall derive the asymptotic formulas of J^n{x) 
(or, Pn{x)) in the exponential region and oscillatory region, respectively. In the following 
lemma, we provide the asymptotic formula in the exponential region of the solution to 
a general class of difference equation. A similar result was obtained by Van Assche and 
Geronimo in [21]. Here, we adopt the approach developed by Wang and Wong in [23] . 



Lemma 2.1. Let 7r„(x) he monic polynomials defined from the following recurrence re- 
lation 

'Kn+iix) = {x- a„)7r„(x) - 6„7r„-i(x) 

with initial conditions vro(x) = 1 and 7ri(x) = x — a^. Here the constants a„ and bn are 
assumed to be polynomials in n and have the following asymptotic behaviors as n oo: 



an = an^ + an 



i^^^ + OinF^^) (2.19) 

6„ = 62^2p ^ ^^2p-l ^ 0(^2p-2^^ ^2.20) 

where a and b > are not identically zero and p > is a positive integer. Let I be 
the smallest convex and closed interval which contains 0, a — 2b and a + 2b, namely, 
I = [a- 2b, a + 2b] if a - 2b < < a + 2b; I = [0, a + 2b] if a- 2b > 0; I = [a - 2b, 0] if 
a + 2b < 0. Rescale the variable x by x = Xn '■= {n + cYy with y & C \ I, then we have 
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the following asymptotic formula for 7r„(x„) 

1/2 



T^niXn) ~ (y) 



{y-a) + ^/{y-ar-Ab^ 
2y 



X exp{n / log[{y — ar^) + \/ {y — arPy — ib'^r'^Pjdr} 
Jo 

, a AlP'r + a{y — ar) 

X exp[ / — =as — / — ■.ar\ 

Jo Viy ~ ^^^y ~ 4&^-5^^ Jo p>y{y - arY - Ab'^r'^ 

/•I 2Br 

X exp[- / — dr]. (2.21) 

Jo Py (y — ary — Ab'^r^[{y — ar) + y (y — ar)^ — Ab'^r'^] 

Proof. Let 

n 

T^nix) = Y[wk{x), (2.22) 

k=l 

we obtain wi (x) = x — Oq and 

Wk+iix)=x-ak ^ (2.23) 

for k > 1. Since x = x„ = (n + o")^?/, we have, for large n 

, . {Xn - Ofc) + {xn - - Abk . , , - 

WkiXn) = ^ (1 + Ek), (2.24) 



where 



= ^^+1-^^- _ + 2ih+i-h) + ixr.-a,)iak^,-a,) 

2v/(^„-a,)2-46,, 2[(x„-afc)2-46,] I J I J 

pa/^P-i ^ ipb'^k'^P-^ + pakP-\ynP - akP) , ^^^_2 



2V(2/r2P - aA;P)2 - 462A;2p 2[(?/nP - aA;P)2 - 462/i;2p] 



+ 0(n-2) (2.26) 



uniformly for = 1,2, ■■ ■ ,n. Here we have made use of the facts that Cfc and bk are 
polynomials in k, and k/n = 0(1) uniformly in /c = 1, 2, ■ ■ ■ , n. On the other hand, 
noting that x„ ~ ynP + S,nP~^ with ^ = pa?/, we obtain 

(xn - afc) + V(^n - - 4fofc = (yn^ + e^^"' - aA;^ - "fc^"') 

+ sJ{ynP + ^raP-i - afcp - akP-^f - Ab^k^P - A[5k^P~^ + ©(n^-^). 
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A simple calculation yields 



where 



1 + 4, (2.27) 



^/{ynP - akpy - Ab^k^P 

- , '^-^^ ^ + 0(n-^) (2.28) 

^J{ynP - akpy - 4b^k^P[{ynP - akP) + ^J{ynP - akpy - Ab^k'^P] 

uniformly for k = 1,2, ■ ■ ■ ,n. By trapezoidal rule, we have the following three asymptotic 
formulas: 

\og[{ynF — ak^) + \J [ynP — akPy — Ab'^k'^P] ~ nplogn + n / log[{y — as^) 
k=i ^0 

+ ^iy-asPy-4b^s^P]ds + log iv - -) + V(^y - (2.29) 

Z Zy 



and 



pasP-^ Apb'^s^P-^ +pasP-\y - asP) 

^^''^ /q 2y/{y-asPy-Ab^s^P ^ 2[{y - asP^ - ib^s^P] 

4b^r + a{y-ar) ^^.30) 



k=l 



2v/(2/-ar)2-46V2 2[{y - ar)^ - Ab^r^] 



and 



y^^fc~ [ — ^ '^^ =ds 
^ Jo ^/{y - asPf - Ab^s^P 



k=l 



2l3s^P^^ 

-as 



y/{y - asPy - 4:b'^s'^P[{y - asP) + ^J{y - asPy - Ab'^s'^P] 

^ fl — [ flfy 

a/(i/ - asPy - Ab'^s'^P Jo P\/{y - ary - Ab^r"^ 

- f = '^-^ ^dr. (2.31) 

Jo Py {y — ary — 4b'^r'^[{y — ar) + y (y — ary — Ah^r"^] 

Combining the above three asymptotic formulas yields the desired formula. □ 
A direct application of the above lemma yields the following result. 
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Proposition 2.3. As n ^ oo, we have 

2n 



-^j exp{nlog(y + + (2n + l)y^F(arcsin(— ), -1)} 

(2.32) 



, \ 1/2 

y + V - 1 
2VF^ 



for complexy bounded away from the interval [—1, 1], where F{ip, k"^) is the elliptic integral 
of the first kind in (12.141) . 

Proof. Choose a = a = 0, 6 = 4, /3 = 0, p = 2, cr=l/2 and replace y by 8y. We obtain 
from flCTD 



J-„(x)~(2n)2"exp[n / \og{y + y^f~t^)dt + _^=rft]( 



/2 



2v/^^ 

(2.33) 



Note that 

/ log(2/ + ^ff^^)dt = t \og{y + Vl/'-t^) + / 
Jo 



2*4 

-dt 



y'-t'{y + Vy'-t') 



log(|/ + VV^) + 2v^ / ^ , ^ 

Jo Vl - + 1 - 

f^Vy 1 

log{y + ^/y^ - 1) -2 + 2^ ^ rfg 

Jo V 1 - 



log(|/ + - 1) - 2 + 2v/^F(arcsin — I - 1), (2.34) 

V 2/ 



and 



, zdt = a/v / ds = A/yF(arcsin , —1). (2.35) 

Combining the above three equations yields (I2.32p . □ 



Once we have the asymptotic formula for J^n{x) outside the interval where the zeros 
are located, using an argument similar to that in the proof of [23| Theorem 3.2], we 
obtain the following result in the oscillatory region. Due to the symmetry of J-'„(x), we 
consider positive x only. 
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Proposition 2.4. Let 6 > be any fixed small number. As n —> oo, we have 

g(n+i)pycos 



J'„(8(n + l/2)'cos^) ~ (2n)'"e 



2n-2n 



sin 2 9 



cos( {n + -){e ^ — Ssm2e(2'4)) 



I / In V COS 6 ^ ,1 1, 
+ sin( {n + -){e 



(2.36) 



for 6 G [5, | — 5], where Bx{a, b) is the incomplete Beta function in (12.151) and the constant 
p is 



p = 2 



1 , , 2v/?r( 



(2.37) 



Finally, we obtain our theorem 

Theorem 2.1. Let u = n + ^. With and U{t) defined m (Q, fl2:T2|) and fl233ll . 

respectively, we have 



r(|) \ /64t/(t)V 



AK.if/(t))f:i§)+A.'(.i^(t))f:^^(^) 



s=0 



uniformly for t in compact subsets of (0, oo). 

Proof. Let ?/ = ''"g^a'' then Sri^y = i^^t = x. Because 

(t/sy. 



fy 1 
n I du ~ (n + -) 



t V64 - s2 V64 - t2 ■ 



we have from (I2.3P and (I2.36P 



Pn{X) 



(2n) 



2n^-2n , /q An+^)p^/t/S 



Kn (64-t2)i 



COS I [n + -) cos 



(n + 



(t/s 



t 764^ 



1, 



+ sm I (n H — ) cos (^^ + - 

2 8 



2' .It V64 - s2 



Recall the asymptotic expansion for Airy function, as x — )■ oo 

Ai(-x) r 



1 ,2 3 TT, 

cos(-x2 - -), 



- -^^sin(^xi-^), 
11 



(2.38) 



(2.39) 



(2.40) 

(2.41) 
(2.42) 



64 \* 1 f2v, 3 TT 



then, we have when t < 8 



Pn{x) = Ci(a:)P„(x) + C2(x)g„(x)~Ci(x)(^^^-^J cos(^y (-[/(t)) 

64 I (2v , , ,,3 TT 

2 



2 



Comparing fl2.40p and (12.431) . we have 



1 /nf) 



This completes the proof of the theorem. □ 

Remark 1.1. Recall the weight function for J-'n(a;) given in (11. 7p . Let a = 0, then we 
have the weight function 

wix) = , X e M. (2.45) 

2[cos(pv/aV2) + cosh(pv^)] 

When a; is large number as in (I2.38p . that is x = z/^t, then we have 
So the main formula in the above theorem can be stated as 



Ai(.it/(t))f;iig)+A.'(.i^(t))f:^^(^) 



« ^^^^ 



(2.47) 

Remark 2.3. In fact, the three-term recurrence relation contains a lot of useful informa- 
tion. It can also give us the asymptotic zero distribution of the rescaled polynomials 
J-'„(8n^x) as n oo. Using the method developed in Kuijlaars and Van Assche [17], one 
can obtain the following limiting measure for J-^„(8n^x) from the recurrence relation (II. 6p 



1 r 1 



TT 



-.ds, xG[-1,1]. (2.48) 
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Since the above theorem is uniformly valid in the neighbourhood of the large transition 
point 8, we have the following asymptotic formula for 

Corollary 2.1. Let i/ = n + |, Kn c^nd w{x) he given in ( 12. 3 p and (12. 45 p . respectively. 



Uniformly for a hounded real numher s, we have, for x 



8v^ 



4 

S Z/3 



w(x)^J^Jx) 



as V 



oo. 



217/6 7]-3/2 ' ' " ■ (2.49) 

Proof. Let t = 8 + 8v^ sv 3 in (I2.47P and recall the asymptotic formula for Uit) in 
(I2.18p . we obtain after some computations 



w(x)2j^Jx) 



2"^/3 + 0(z/^i)) Ai(s) + 0(z/-i) 



(2.50) 
□ 



4^2 713/2 

This above formula proves our corollary. 

Remark 2.4. We may rewrite the polynomial on the left hand side of (12.491) into an or- 
thonormal one. Note that if Vn{x) satisfies a three-term recurrence relation in (II. 8p . then 
Vn{x) := {]Xk=iPk)~'^Vn{x) is the corresponding orthonormal polynomial and satisfies 
the following recurrence relation 

xVn{x) = ^/K^lVn+l{x) + anVn{x) + ^/KVn^l{x). (2.51) 

Let J^n{x) be the orthonormal Chen-Ismail polynomials, then (I2.49P can be rewritten as 



w{x)^Fn{x) 



5 

V 6 



27/3 _ 



as z/ — > oo. 



(2.52) 



Ai(s) + 0(i^~3 

From the above corollary, an asymptotic approximation of extreme zeros of J^n{x) can 
be easily obtained by using the zeros of the Airy function Ai(x). To get it, we need the 
following result of Hethcote [12] . 

Lemma 2.2. In the interval [a — p, a + p], suppose that f{t) = g{t) + e{t), where f{t) is 
continuous, g{t) is differentiate, g{a) = 0, m = min \g'{t)\ > and 

E{t) = max \e{t)\ < min \g{a - p)\, \g{a + p)\. (2.53) 

Then, there exists a zero c of f(t) in the interval such that \c — a\ < ^. 

We have the following approximation. 

Corollary 2.2. Let Xn,k he zeros of J^n{x) such that Xn,i > Xn,2 > ■ ■ ■ > Xn,n, and ak he 
zeros of the Airy function Ai(x) in the descending order. Then we have for fixed k and 
large n 



Xn,k = 8u^ + 8\/2 Ofcl 



,4/3 



0{u 



2/3N 



(2.54) 



where u = n + ^ . 



Proof. A combination of Corollary 12.11 and Lemma [2.21 immediately gives us the result. 

□ 
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3 Berg— Valent polynomials 



Now we are going to study some birth and death process polynomials Qnix). Before we 
derive the Plancherel-Rotach asymptotics for Berg- Valent polynomials Qnix) with the 
rates in (11.31) . we can use the chain sequence method to get the bounds for the largest 
and smallest zeros as we did in Proposition 12.11 

Proposition 3.1. Let Xn,k be zeros of Qn{x) such that Xn,i > Xn,2 > ■ ■ ■ > Xn,n, then we 
have the following hounds for all n > 1 

Xn,i < 2^V - 2^V + 35 ■ 2^n'^ and Xn,n > 4.29. (3.1) 

Proof. Recall the recurrence relation of Qn{x) in (II. 3p and choose c„ = 1/4 in Theo- 
rem 11.11 Note that the solutions Xn and i/n in (11.111) are monotonically increasing and 
decreasing with n, respectively. Then the result follows. □ 

3.1 Difference equation method 

Like what we have done in the previous section, we introduce Qn{x) = (— l)"-ft'„p„(x) 
with 

:= (3.2) 

r(^)2r(^) 

to arrive at the standard form. Note that 

(4n-l)(4n)2 



Kn-i (4n + 2)2(4n + 3)' 
then the recurrence relation (II. ip with (II. 3p reduces to the following standard form 

P„+l(x) - {AnX + Bn) Pn{x) + Pn-l{x) = 0. (3.3) 

As n — oo, the recurrence coefficients An and -B„ satisfy the following expansions 

An = T TP n > — , Bn = 7 > — 3.4 



with e = 4, 



and 



an = , ai = (3.5) 

256 256 ^ ^ 



/3o = -2, (3i=(32 = 0. (3.6) 
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The characteristic equation for (13.31) is 

- {aot + f3o)X + l 







(3.7) 



with ao and f3o given in (13. 5p and (13. 6p . respectively. The roots of this equation coincide 
when aot± + /3q = ±2, which gives us two transition points 



,10 



0. 



(3.8) 



Near the large transition point t+, we get the Airy-type asymptotic expansion as in the 
previous section. But since the small transition point t_ is located at the origin, we 
obtain the Bessel-type expansion in its neighbourhood. For an explanation why the case 
t_ = is so special to give us the Bessel type asymptotics, one may refer to discussions 
in [6]. Then we have the following two different types of expansions. 



Proposition 3.2. When n is large, Pn{x) in (13.31) can be expressed as 

Pn{x) = Ci{x)Pn{x) + C2{x)Qn{x), 



(3.9) 



where Ci{x) andC2{x) are two n-independent functions, P„(x) andQn{x) are two linearly 
independent solutions of (13.31) satisfying the following Airy-type asymptotic expansions 
in the neighbourhood oft^ = 2^^ 



A,(.ic/w)f;^+A,vic/w)E'^'"'' 



s=0 



s=0 



(3.10) 



and 



Q^iv't) = 2^/2 



Uit) 



t{t 



Bi{ulu{t))Y, 



As{U) 



s=0 



V 6 



B^{ulu{t))Y, 



s=0 



Bs{U) 



Here p = n + \, the leading coefficients are given by 



Bo{U)=0 



and U (t) is defined as 



Imt)] 

l[-u{t)] 



y/2t^ 



1 1, 



cos 



-1 



29 



log 



t-2^ + ^t{t - 210) 



29 

^1 * (-,-) 



(3.11) 
(3.12) 

t > 2^° (3.13) 
< t < 2^°; (3.14) 



see (I2.15P for the definitions of Bx{a,b). 



15 



Proof. With 

ait+ + /3i _ 1 
' (2 - f3o)9 ~ 4 

our results follow from the main Theorem in 



To 



and u = n + To, 



(3.15) 
□ 



Proposition 3.3. Similarly, when n is large, Pn{x) in (13 .Sp can be expressed as 

p4x) = C*(x)(-1)"P:(x) + C;{x){-irQ*n{x) (3.16) 

where C*(a;) and C2{x) are two n-independent functions, Pn{x) and Q*^{x) are two lin- 
early independent solutions of (13. 3p satisfying the following Bessel-type asymptotic ex- 
pansions in the neighbourhood oft^ = 



and 



U*{t) 
t(2io - 1) 



U*{t) 
t(2io - 1) 



B*{U* 



s=0 



(3.17) 



W.{uU*Ht)) J2 + W-^{vU*Ht)) J2 



Here u = n + ^, 

the leading coefficients are given by 

i*(f/) = l, B;iU) = Q 

and U*{t) is defined as 

itls)\ 



(3.18) 
(3.19) 

(3.20) 



, , 29 - t + Jtit - 210) 
ds - log — ^TT^ t < (3.21) 



o t (-,-) — cos 

8 2^ 



29 



29 



< t < 2^"; (3.22) 



see f l2.15p /or i/ie definitions of Bx{a,b). 



Proof. With the asymptotic expansions for and i?„ in (13.40 . our results follow from 
the main Theorem in [B]. □ 



Remark 3.1. Note that the integral in (I3.2ip can be written as a hypergeometric function 
as follows: 



't \/S{S 



210) 



2 2^iM'2'4'' 
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The two functions U(t) and U*(t) are both monotonically increasing function at 2^° and 
0, respectively. In face, from their definitions in the above two propositions, we have the 
following asymptotic formulas 

f - 2^° 

Uit) = -^ + 0{t-2'y, ast^2i°. (3.24) 

and 

U*{t) = ^ + 0{f) ast->0. (3.25) 
3.2 Determination of C\{x) and C^ix) 

For convenience, we state a special case of Lemma 12.11 which will be used in determining 
C\[x) and C2(x) for all three types of birth and death polynomials considered in this 
paper. 

Lemma 3.1. Assume 7rn{x) satisfies the following recurrence relation 

Tln+lix) = [X- (A„ + /i„)]7r„,(x) - Xn-lflnTln-lix) (3.26) 

with 7ro(x) = 1 and vri(x) = x — {Xq + Hq). Here, A„ and fi^. are assumed to he polynomials 
in n and satisfy the following asymptotic formulas as n ^ oo: 

A,, = b{nP + unP-') + 0{nP-^) (3.27) 
= h{nF + vn^-^) + 0{nF-^), (3.28) 

where b > andp > 0. Rescale the variable x by x = Xn '■= {n + aYy with y E C \ [0,46], 
then we obtain the following asymptotic formula 



vr.(x.)~(^r(l-4Vy)-^/' 



1 + y/l-Ab/y 



X expj / — (»- + o')p 

(3.29) 



Proof. We will show that fl3.29p is a special case of f l2.2ip with a = 2b, a = b{u + v) 
and /3 = b'^{u + v — p). Due to the fact that a = 2b, we have following formulas for the 
integrals in fl2.2ip 



dr = I — dr = — . 



2 A/(y — ar)^ — Ab'^r'^ Jo 2 \/y'^ — 2ayr 
and 

'1 46^r + a{y — ar) ay r 1 i V 



dr = — m rdr = - lo£ 



2[{y - ary - Ab'^r^] Jq 2[y^ - 2ayr] 4 y - 2a 
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and 



and 



a 



Pa/ [y — ary — Ab'^r'^ 



dr 



ail - y/l - 2a/y) 
pa 



2/3r 



p^J {y — arY — 4:b'^r'^[{y — ar) + \/ {y — arY — ^h'^r'^] 
4/3 



dr 



■ [log 

pa?' 1 + - 2a/y 



1 - - 2a/|/, 



Therefore, we have from f l2.2ip 



^nixn) ~ (y)"(l - 2a/yr'/' 



1 + 2a/y 



-1 1- 



4l_ 



X exp{n / log[{y — as^) + \/ y'^ — 2aysP]ds + / 
io Jo 



pay 



ds} 



X exp[(^ - - + ^)(1 - A/l-2a/2/)]- 
2 pa pa^ 



(3.30) 



Furthermore, since 



(y - asP) + ^/y^ - 2aysP = (Vv + ^y-^as^fj^, 



we have 



I log[(?/ - as^) + ^/y^^^2aysP]ds = - log2 + 2 /" logf^/y + - 2as?']rfs. 
Jo Jo 



According to integration by parts, the right-hand side becomes 

-1 ^-Vy- 2asP 
Vy - 2asP 

ds. 



log2 + 21og[v/^+ Vy-2a] +p 



ds 



log — p + 



^ p^/y 



Therefore, we obtain 



1 + ^1- 2a/y 



- 2asP 



2n+l+- 



^ (n + (t)p 



,1 a 2(5 



ds + {-- — + ^)(1 - - 2a/y)}. (3.31) 



v^l - 2asP/y 2 pa pa^ 
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Note that 



1 a 2/3 1 u + v ^u + v — p ^ 



2 pa pa^ 2 2p 
The asymptotic formula fl3.29p for 7r„(x„) follows 



2p 



□ 



Since the birth and death process polynomials Qn{x) defined in (11.11) is related to the 
monic one as follows 



Qnix) 



n-1 



k=<d 



-1 



(3.32) 



we have the following results for Berg-Valent polynomials from the above lemma. 
Proposition 3.4. Let x = Xn '■= {n + 1/4)H with t G C \ [0, 2^°] , we have as n ^ oo, 

2n+l/2 



p^{x^) ~ 2-«'^-5/2(n + 1)^(1 - 2^7^)-!/^ 



1 + V^l - 210/^ 



X exp{^^^jy^F(arcsin(2^7t)V4| _ 1)}. 
Proof. Recall the definition of Kn in (13. 2 p , we have 



(3.33) 



r(^)r(^)2 



\n. 2 



1 ' 

4 



which yields Qn(a^) = {—^Y KnPn{x) ~ (— 1)"— -TPn(3^)- Moreover, by applying Stirling's 
formula, we obtain 



fc=0 



28"r(l/4 + n)r(l/2 + n)2r(3/4 + n) 

r(i/4)r(i/2)2r(3/4) 



It is readily seen that 



r(i/4)r(i/2)2r(3/4) 

28"(27r)2(n/e)^" 

r(i/4)r(i/2)2r(3/4) 

2«'^2y2(n/e)^". 



n-l 



7r„(x) = {-irT\\,Qn{x) ~ 2«"+=^/2(n/e)^" jp^{x) 



A:=0 



Since '/r„(x) satisfies the asymptotic formula in (I3.29P with a = 1/4, p = A, b = 2^,u = 2 
and w = 0, our result follows from the above formula. □ 
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Again, using an argument similar to that in the proof of [231 Theorem 3.2], we obtain 
the following result in the oscillatory region. 

Proposition 3.5. Let x = Xn '■= {n + 1/4:)H with t = 2^^ cos'^ 6,6 E [<5, f — S], we have 
as n CO, 



1 1 

Pn{xn) ~ V2{n + -)-——— J exp 
(/ t) i 



(n 



TT 1 

X cos{- - (2n + 1/2)6 + (n + -) 



{t/s 



1 1/4 



V2 



10 



s — 



-.ds} 



where 



P 



du 



1 Ui\/u^ — u 



r(|) 



(3.34) 
(3.35) 



Then we get our theorem in the interval containing t+ = 2^^. 

Theorem 3.1. Let u = n + I. With and U{t) defined m (O, fl3:T3|) and ^1^ . 

respectively, we have 



-l)"A':„V27r vti exp 



^ / v^r(i) A ( U{t) 
4 exn — T- — -^vti — 



\^25/2r(| 



X 



Ai(.if/(t))^^ + Ai'(.if/(t))5^ 



t{t - 210) 
Bs{U) 



s=0 



(3.36) 



uniformly for t in compact subsets of (0, oo). 

Proof. When t < 2^°, recall the asymptotic expansion for Airy function in fl2.4ip and 
(^M> again, we have from f lXTU]) and fIXTT]) 



Pnix) = Ci(x)P„(x) + C2{x)Qnix) ~ Ci(x) 
29/2 1 



29/2 



r^cos — (-f/(t))2 - - 

[t(2i0-t)]i V3^ 4 



2i/ 



3 TT 



- Co(x) 



CAx) 



2z^ / , - - 3 vr 
, — -f/ t 2 _ _ 

[t(2io_t)]i V3^ 4 



sm 



29/2 



1 



[t(2i0-t)]3 ^/7r 



COS I z/(cos ^ 



t -2 



C2(X) 



29/2 



[t(2i0-t)]3 



sm z/ cos 



-1 



29 
t-29 



29 



M=d.) - ^ 

s(2io - s) 4 



10 



Comparing the above formula and fl3.34p . we have 



Ci(x) 



24 



exp X * 



iv^r(i; 



25/2 r(|) p 



Co(x) = 0. 



This completes the proof of the theorem. 



(3.37) 
□ 
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Remark 3.2. Like in Remark 12.31 we can obtain the following limiting zero distribution 
for Q„(2^°n^x) from the recurrence relation (11. 3p 



^/x{s^ - 



X) 



-.ds, X G [0, 1]. 



(3.38) 



Since the above theorem is uniformly valid in the neighbourhood of the large transition 
point 2^°, we have the following asymptotic formula for Q„(x). 

Corollary 3.1. Let u = n + j and Kn he given in (13. 2p . Uniformly for a hounded real 
numher s, we have, for x = 2^^iy^ + 2^°v^s u~ 



X 4 exp 



25/2 r(| 



^' xi 1 QJx) 



16^ 



ki{s) + 0{u 



as V ^ oo. 

(3.39) 



Proof. Let t = 2^° + 2^^^/l sv 3 in (I3.36P and recall the asymptotic formula of U (t) in 
(I3.24p . we have after some computations 



X 4 exp 



v^r(i) 
25/2 r(|) 



X4 QJx) 



32\/2 



Z/6 



This above formula proves our corollary. 



Ai(s) + 0(z/~3) . (3.40) 



□ 



Remark 3.3. We can put the left-hand side of (13.390 into its orthonormal form as in 
Remark 12.41 Let Qnix) be the orthonormal Berg-Valent polynomials, then (I3.39P can be 
rewritten as 



X 4 exp 



25/2 r(|; 



X* Qnix) 



_ 11 
U 6 

16v^ 



Ai(s) + 0(i/"2 



as I' oo. 



(3.41) 



Unlike (12. 52 p . we don't have w{x) on the left hand side of the above equation because 
the weight function w{x) is unknown for the Berg-Valent polynomials. But it seems 
reasonable to conjecture that at least one of the weight functions for the Berg-Valent 
polynomials should behave like 



X 2 exp 



v^r(i) 1 



23/2 r(|) 



X* 



as X — >■ oo. 



(3.42) 



We shall revisit this issue in §6. 

We have the following approximation. 
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Corollary 3.2. Let Xn,k be zeros of Qn{x) such that Xn,i > Xn,2 > ■ ■ ■ > Xn,n, and be 
zeros of the Airy function Ai(x) in the descending order. Then we have for fixed k and 
large n 

Xn,k = 2^V' + aku'""^ + 0(z/S/3), (3.43) 

where u = n + ^. 

Proof. A combination of Corollary 13.11 and Lemma [2.21 immediately gives us the result. 

□ 

We also get our theorem in the interval containing t- = 0. 

Theorem 3.2. Let u = n + I. With and U*{t) defined m flX^ and flX^ . 

respectively, we have 



Qn{iy^t) = K^VtI Z/2t4 exp 



1,25/2 r(| 



-ut* 



U*{t) 
t(2io - t) 



X < sm 



cos 



V25/2r(f 



-Vti 



s=0 s=0 



(3.44) 



uniformly for —oo<t<M <2^°. 

Proof. When < t < M, we have from (13X7]) and fl3A8|) 

p^{x) = i-irci{x)p:ix) + i-irc;{x)Q:ix) 



-1)"29/Vi 



U*{t) 



CI (x) ji {uu* Ht)) + c; {x)Wi{i^u*^ (t) ) 



_t(2iO-t), 

As Pn{x) is real, we may choose Cl{x) = Ci{x) + iC2{x) and 6*2(2;) = C2{x) such that 

U*{t) 



p„(x) ~ (-1)"29/V^ 

nr)5 ' 



-1)"2 



t(2io - 1) 

1 

^ 1 



Ci (x) Ji (i^f/* 2 (t) ) + C2 (a;)n (i^t/* 5 it) ) 



Coix) cos 



t(2i0-t); V? 

(t/s)i 



C*i(x) sin 



(t/.)3 _i/29-t 
=as — cos 



10 



29 



Vs(2io - s) 



as — cos 



29 



(3.45) 



Since 



210 



V2i0s - s' 



-.ds = 2-'/^pt-^- f 
' Jo 



(t/sy/^ 



V2 



10 



-ds 



s — s^ 
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and 20 = -n — cos'^ 



-1 I 2'>~t 
29 



we rewrite f l3.34p as 



Pn(x) ~ V2{n + ^) : 

r ^ 

X COS-^ \- UTl — UCOS 

^ 4 



exp 



-1 



(n 



)9 



V-2 "''-^".exp 

(2iot-t2)i 



X COS 



29 



,1/4 



V2i0s - s' 



-.ds} 







v/s(2io - s 

Comparing f l3.45p and the above formula, we have 



29 



2-5/2 



Clix) 



C2{X) 



TT X4 

2"97^ 

— j_ 

TT X4 



exp X* 



4 Ar(i; 



29/2 



exp I X 4 



25/2 r(|) 



vrr(i) 



sm x4 



25/2 r(|) 



cos X4 



v^r(i) 
25/2 r(f) 

25/2 r(|) 



This completes the proof of the theorem. 



(3.46) 



(3.47) 



(3.48) 



□ 



Remark 3.4. To get the asymptotic formula for t < from (13.441) . one need to consider 
the value t±ie and take the limit as e — )■ 0. The interested readers may compare the for- 
mula (13.441) with the Bessel-type asymptotic expansion for the Laguerre-type orthogonal 
polynomials in [221 Eq. (2.15)]. 



4 Conrad— Flajolet polynomials I 

We can easily get a bound for the largest zero for Qn{x) from Theorem II. 1[ 

Proposition 4.1. Let Xn,k be zeros of Qn{x) such that Xn^i > Xn,2 > ■ ■ ■ > Xn,n, then we 
have the following hounds for all n > 1 

Xn,i < lOSn^ + lOScn^ (4.1) 

Proof. Recall the recurrence relation of Qn{x) in (11.41) and choose c„ = 1/4 in Theorem 
11.11 Then the result follows. □ 

Note that we don't provide a bound for the smallest zero in this case. The reason is 
that, although the chain sequence method can give us an estimation, it is not accurate. 
In fact, we know that all the zeros of Conrad-Flajolet polynomials should be positive, 
but Theorem II. II can only give us > — 6n — 3. For the similar reason, we only consider 
the largest zero in Proposition 15. H too. 
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4.1 Difference equation method 

As usually, we introduce Q„,(x) = {—l)^KnPn{x) with 



- "+§ + l ^2 



(4.2) 



to arrive at the standard form. Note that 

Kn+i (3n + c)2 



Kn-i (3n + c + 2)2' 
then the recurrence relation (11. ip with (ll.4p reduces to the following standard form 

Pn+lix) - (AnX + Bn) P„,(x) + Pn^l{x) = 0. (4.3) 

As n — > oo, the recurrence coefficients satisfy the following expansions 



with ^ = 3, 



and 



n 



n+l 



1 

" 27' 



(An + Atn)-ft'n 



n+l 



s=0 



ai 



c+1 
27 



/3o 



/3i = 0, 



/32 



9' 



(4.4) 



(4.5) 



(4.6) 



Then the roots of the characteristic equation (13.71) coincide when aot± + /3o = =t2, with 
ao and /3o given above. This gives us the following two transition points 



108, 



0. 



(4.7) 



These transition points are similar to what we have for Berg-Valent polynomials. So the 
following two Airy-type and Bessel-type expansions follow. 



Proposition 4.2. When n is large, Pn{x) in (14. 3 p can be expressed as 

Pn{x) = Ci{x)Pn{x) + C2{x)Qn{x), 



(4.8) 



where Ci{x) andC2{x) are two n-independent functions, Pn{x) andQn{x) are two linearly 
independent solutions of (14.31) satisfying the following Airy-type asymptotic expansions 
in the neighbourhood of t^ = 108 



m 

t{t- 108) 



Ai{ulu{t))J2 



V 6 



k^{vlu{t))Y^ 



Bs{U) 



V ' 6 



(4.9) 



24 



and 



Qn{uH) = 3^6 



Uit) 



tit - 108) 



Bi(.ic/w)f:Mj)+Biv!c/w)f:^'(''> 



s=0 



V 6 



i/ere z^ = n + i/ie leading coefficients are given by 



Ao{U) = l, Bo{U)=0 



and U (t) is defined as 



l[-u{t)]l 



y2i ^ ,11, , t - 54 + .yt{t - 108) 

Oi 108 -,- — log 

6 ^- — ^2' 3^ ^ 54 



cos 



t -54 
54 



-r^^-w^2'6^ 



(4.10) 
(4.11) 

t > 108, (4.12) 
< t < 108; (4.13) 



see 02.151] for the definitions of Bx{a,b) . 
Proof. With 



To 



+ /3i c + 1 



(2 - /3o)^ 3 
our results follow from the main Theorem in |2^. 



and u = n + To, 



(4.14) 

_ □ 

Proposition 4.3. Similarly, when n is large, Pn{x) in f l4.3p can 6e expressed as 

p„(x) = Cr(x)(-l)"P:(x) + C2*(x)(-l)"g:(x) (4.15) 

where Cl{x) and C2{x) are two n-independent functions, P^ix) and Q^{x) are two lin- 
early independent solutions of fl4.3p satisfying the following Bessel-type asymptotic ex- 
pansions in the neighbourhood oft_ = 



Pn{yH) = 3^6 
and 



U*{t) 



t(108-t) 



U*{t) 
t(108-t) 



s=0 



s=0 



A:iu*] 



(4.16) 



s=0 



s=0 



Here u = n + Wci{x) is defined in (13.191) . the leading coefficients are given by 



(4.17)' 



i*(t/) = l, 5*(f/) = 



(4.18) 



25 



and U*(t) is defined as 



° , , 54 - 1 + Jt(t - 108) , 
^' ' ds -loK ^ t<0 (4.19) 



6 ^'6 2 
see f l2.15p for the definitions of Bx{a,b) 



v/s(s- 108 

„ ,11, _i /54-t 
i3_(_(-, -) — cos 



54 



54 



< t < 108; (4.20) 



Proof. With the asymptotic expansions for and i?„ in (14.40 . our results follow from 
the main Theorem in [B]. □ 



Remark 4.1. Again, the integral in f l4.19p can be written as a hypergeometric function. 

(4.21) 



° (t/s)^ , 3/-, ,1 ^ ,1 1 7 , 
^s(s - 108) 6 2 6 

U{t) and f/*(t) are two monotonically increasing function at 108 and 0, respectively, with 
the following asymptotic formulas 



[/(t) = ^-^S + 0(t - 108)2, as 108. 



54^ 



and 



At 

U*{t) = — + 0{t^) ast^O. 



(4.22) 



(4.23) 



4.2 Determination of CAx) and C2(x) 



According to Lemma [3.11 we have the following result. 

Proposition 4.4. Let x = Xn '■= {n + (c + l)/3)'^t with t G C \ [0, 108], we have as 
n — )■ oo, 



, , r((c+l)/3)r((c + 2)/3)2t" ^^ in«m-i/4 



33n^c+l(2vr/n)3/222/3 
^ 3n + C + 1 



1 + a/1 - 108/t 



2n+l+ 



(2c-l) 



X exp{ / — 
Jo \/l 



108sVt 



(is}. 



(4.24) 
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Proof. Recall the rates given in (11 ■4p and Kn defined in (14. 2p . By Stirling's formula, we 
have Kr, ~ 2^/^?i"^/^ and 



fc=0 



33"r((c + l)/3 + n)T{{c + 2)/3 + nY 
r((c+l)/3)r((c + 2)/3)2 

r((c+i)/3)r((c + 2)/3)2 

33n^c+5/3 ^2vr /n) (n/e) 
r((c+l)/3)r((c + 2)/3)2- 



It thus follows from (I3.32p that 



n-l 



k=0 



33n^c+i (^2vr /n)3/2 (n /e)3"22/3 
r((c+l)/3)r((c + 2)/3)2 



Since 7r„(a;) satisfies the asymptotic formula (13.291) with p = 3, b = 27, u = c + 5/3 and 
V = c + 1/3, our result follows from the above formula. □ 

Again, using an argument similar to that in the proof of [23l Theorem 3.2], we obtain 
the following result in the oscillatory region. 

Proposition 4.5. Let x = x„ := (n + (c+ 1)/3)H with t = 108cos^6',6' G [5, f - 6], we 

have as n ^ CO, 



PniXr, 



3^+ir(£±i)r2(£±2^ 



27/67r3/2(108t-t2)l/4 (^3^)2^ 



— exp {n 



,Ti , (2c + 2) 
xcos{--(2n+ ^ 3 



c+ L 



3 ^408^ 

l°«_(t/£)V^ 



where 



P 



du 



1 U i — u 



v^r(|: 



:c/s}, (4.25) 



(4.26) 



Then we get the Airy-type expansion in the interval containing t+ = 108. 

Theorem 4.1. Letv = n + ^. With Kn and U{t) defined in IKW\ and flCTj) . 

respectively, we have 



3C+I -p ^ c+l ^ p2 i' c+2 ^ 



v^r(i; 



3 ' 3 ' 

n 7 2c-i 2c-i exp . -jAT-n/c; 



X 



Ai(.if/(t))5^^ + Ai'(.if/(t))5^ 



ut^ 

Bs{U) 



U{t) 



t{t - 108) 



U 6 



(4.27) 



uniformly for < 6 < t < 00. 
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Proof. When t < 108, recall the asymptotic expansion for Airy function in f l2.4ip and 
fl2.42p again, we have from ( I4l9|) and fl4.10p 

3^6 



Pn{x) = Ci{x)Pn{x) + C2{x)Qn{x) ~ Ci{x) 



[t(108-t)] 



4 x/TT 



2u 



IT 



COS 



C2ix 



3^6 



, ^sin( — (-[/(t))t - - 
[t(108-t)]i V3^ 4 



2i^ 



3\/6 



1 



C2(X) 



[t{108-t)]i 
3^6 



cosi z/(cos ^ 



[t(108 -t)]3 \/7r' 



sm 1/ cos 



t-54 
54 

t - 54 
54 



t Vs(108-s 
t Vs(108 - s) 4 ' 



Comparing the above formula and fl4.25p . we have 
_ 3-V^r((c+l)/3)r^((c + 2)/3) 



This completes the proof of the theorem. 



108' 



r(i; 



C2(a;) = 0. 



(4.28) 
□ 



Remark 4.2. Like in Remark 12.31 we can obtain the following limiting zero distribution 
for Qn{108n^x) from the recurrence relation (11. 4p 



--ds, X e [0, 1]. 



, . (4.29) 

^ y x(s^ — x) 

Since the above theorem is uniformly valid in the neighbourhood of the large transition 
point 108, we have the following asymptotic formula for Q„(x). 

Corollary 4.1. Let u = n + ^ and Kn he given in (14. 2p . Uniformly for a hounded real 
numher s, we have, for x = lOSu^ + 54v^r8 s 



2c-l 

X 6 exp 



v^r(|; 
3^r(| 



-X3 Qnix) 



3c-ir(£±i)r2(£±2)j,i 



Air 



Ai(s) + 0(z/"3 



as u oo. 



(4.30) 



Proof. Let t = 108 + 5A\/T8si' 3 in (I4.27P and recall the asymptotic formula for U{t) in 
(I4.22p . we have after some computations 



2c-i / v^r(i) 



X 6 exp 



XZ/6 



3^r(|) 

Ai(s) + 0(z/ 



X3 Qn{x) = {-irKn 



„^ 3^+^r(^)r2(^) / 1 



267r 



3326 



(4.31) 



This above formula proves our corollary. 



□ 
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Remark 4.3. Again we can put the left hand side of f l4.30p into its orthonormal form. Let 
Qnix) be the orthonormal Conrad-Flajolet polynomials I, then fl4.30p can be rewritten 
as 



2e-l 

X 6 exp 



3^r(|) 



x-^ QJx) 



21 35 



Ai(s) + 0(i^"3 

as oo. (4.32) 



According to the above formula, again it seems reasonable to conjecture that at least one 
of the weight functions for the Conrad-Flajolet polynomials I should behave like 



2^ / v^^r(i) 



X 3 exp 



3r(i) 



-Xi 



as x — )■ oo. 



(4.33) 



We have the following approximation. 



Corollary 4.2. Let Xn,k be zeros of Qn{x) such that Xn,i > Xn,2 > ■ ■ ■ > Xn,n, and be 
zeros of the Airy function Ai(x) in the descending order. Then we have for fixed k and 
large n 

Xn,k = 108i/^ + 54v^ OfcZ/^/^ + 0(i/^/2), (4.34) 

where u = n + ^ . 

Proof. A combination of Corollary 14. II and Lemma [2.21 gives us the result. □ 
We also get the Bessel-type expansion in the interval containing t_ = 0. 

Theorem 4.2. Letu = n + ^. With and U*{t) defined m (g^]), flCTj) and flCTj) . 

respectively, we have 

'ic+i-r(c±i\y2/c+2\ i-c 



r ;o 2c-l 

25/% t— 



exp 



X < sm 



v^r(|) i c 

^5 — — X'i TT 

2iv^r(|) 3 



f v^r(|: 

— 1= — ^^t^ 

V3^r(|) 



U*{t) 



t(108-t) 



s=0 



— COS — — — ^—X3 TT 

2§V3r(|) 3 



Y.iuU*Ht))f2 



Al{U* 



(4.35) 



uniformly for —oo < t < M < 108. 
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Proof. When < t < M, we have from (Hl6|) and (HTTD 

pM = {-irci{x)p:{x) + {-irc;{x)QUx) 



-l)"3y6i^5 



U*{t) 



CI (x) Ji {uU* 2 (t) ) + C; (x) Wi {vU*-^{t)) 



t(108-t) 

Note the following asymptotic expansion for Bessel function, as x — ?■ cxd 



f 2 



f 2 



■^^^"^^^^^^J cos(x----), y^{^)-[-j sin(x----). (4.36) 



As Pn{.x) is real, we choose Cl{x) = Ci{x) + iC2{x) and (x) = C'2(x) such that 

1 



_t(108-t) 
,6^3/ 1 



C'i(x)Ji (z/f/*2(t)) + C'2(a;)Fi(z/f/*2(t)) 



X <^ Ci{x) cos z/ 
C*2(x) sin 



n \t{108-t) 



V 



v/s(108-s; 



:(is — cos'^ 



1 



54 



12 



-TT 



v^s(108-s) 



— cos'^ 



1 /^54-t 
54 



12 



-TT 



(4.37) 



Since 



108 



t Vl08s - s2 2iy3r( 



2 /::r„,5, 
6> 



Vl08s - s2 



and 26' = TT — cos ^ (^|4^) we rewrite fl4.25p as 

Vn{x) 



3c+ir(£±l)r 2(£±g) 1 

27/67r3/2(108t 



r vr _i /54-t 

X COS-^ h Z/TT — Z/ COS 

^ ^ V 54 



2iV3r(|) io Vl08s - s2 



ds} 



3^+ir(£±i)r2(£±2 



2^6,^3/2 x^(l08t-t2)^'^n408^'^'^ 



X 



X COS 



z/ 



(t/s)3 , , /54-t 

as — cos 



Vs(108-s) 



54 



-TT 



v^r(i; 



,c 1, 



12 23v^r(|) ^3 2^ 

(4.38) 
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Comparing fl4.37p and the above formula, we have 

Ci(x) = ^ 2c 1 exp x3 sin ^ x3 - -vr , 4.39 



3-v^r(^)r^(^) / i v^^rxl) \ / v^^r(|) , c 

— , , ^ — exp x3 — — — ^ cos — — x3 



This completes the proof of the theorem. □ 

5 Conrad— Flajolet polynomials II 

Again we use Theorem 11.11 to get a bound for the largest zeros for Qn{x) as follows. 

Proposition 5.1. Let Xn,k be zeros of Qn{x) such that Xn,i > Xn,2 > ■ ■ ■ > Xn,n, then we 
have the following bounds for all n > 1 

Xn,i < lOSn^ + lOScn^ (5.1) 

Proof. Recall the recurrence relation of Qn{x) in (11.51) and choose c„ = 1/4 in Theorem 
11.11 Then the result follows. □ 

5.1 Difference equation method 

As what we have done before, we introduce Q„,(x) = (— l)"-ft'„p„(x) with 

n+l+l n+f+1 .2 
K — ^ 2 Ml 2 ) 

to arrive at the standard form. Note that 



Kn+i {3n + c - l){3n + c 



,2 



Kn-i (3ra + c+l)2(3n + c + 2)^ 



then the recurrence relation 01. ip with ( II. 5p reduces to the standard form ( 14. 3p . As 
n — 7- oo, the recurrence coefficients satisfy the following expansions 



with 6 = 3, 



1 2c+l , ^, 

«o = ^, «i = (5.4) 
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and 



/3o 



/3i=0, 



36' 



(5.5) 



Since ao and /3o are the same as those in (14. 5p and (14. 6p . we get the same transition 
points as follows 

t+ = 108, t_ = 0. (5.6) 

Again we get two different types of asymptotic expansions near and t_. Due to 
the similarity between the current case and that in Sec. HI we get the same Airy-type 
expansion as that in Proposition 14.21 except that the constant u = n + in this case. 

Because the higher order coefficients in (15. 4p and (15. 5p are different with those in 
(14. 5 p and (14. 6p . we have a corresponding difference in the order of Bessel functions. The 
Bessel-type expansion near t_ is given as follows. 



Proposition 5.2. When n is large, Pnix) in (14. 3 p can be expressed as 

Pn{x) = ci{x){-irp:{x) + ci{x){-irQi{x) 



(5.7) 



where Cl{x) and C2{x) are two n-independent functions, Pn{x) and Qni^) are two lin- 
early independent solutions of (14. 3 p satisfying the following Bessel-type asymptotic ex- 
pansions in the neighbourhood oft^ = 



P*{vh) = 3^6 i/^ 
and 

Ql{vH) = u'^ 



U*{t) 



t{108-t) 



U*{t) 
t(108-t) 



s=0 



s=0 



W.{uW-^{t))Y, 



Al{U* 



W.{vW-^{t))Y, 



(5. 



s=0 



s=0 



Here u = n + Wa{x) is defined in (I3.19p . the leading coefficients are given by 



(5.9)- 



i*(t/) = l, 5*(f/) = 0. 
and U*{t) is defined the same as (14.191) and (14.201) . 



(5.10) 



Proof. With the asymptotic expansions for and i?„ in (15.30 . our results follow from 
the main Theorem in [B]. □ 
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5.2 Determination of Ci{x) and C2{x) 
According to Lemma [3. 11 we have the following result. 

Proposition 5.3. Let x = x„ := (n + (2c + l)/6)^t with t G C \ [0, 108], we have as 
n -> oo, 



r((c+ l)/3)2r((c + 2)/3)t", ,,-1/4 
Vn[x^) - \„JJ,,l^_)^^,JJJ (1 - 108/t)-^/^ 



X exp{ 



33n,^c+l/2(27r/n)3/225/6 

^ 3n + c + 1 /2 



'I - 108/t 



2n+l+ 



(2c-2) 



1 - 108s3/t 



ds}. 



(5.ii; 



Proof. Recall the rates in (11. 5p and Kn defined in (15. 2p . By Stirling's formula, we have 
~ 2^/^-5/6 and 



n-l 
A;=0 



33"r((c + l)/3 + n)2r((c + 2)/3 + 

r((c + i)/3)2r((c + 2)/3) 
r((c+i)/3)2r((c + 2)/3) 

33n^c+4/3 (27r/n) ^'"^ {u / cf"^ 

r((c+i)/3)2r((c + 2)/3)- 



Then from (I3.32p . we have 



n-l 



vr„, X 



TT^ .X 33"n'^+i/2(27r/n)3/2(n/e)3"25/6 ^ , 

Kn II AfcP„(x) T.//„ , ^^ ,n\'>T^n. , ON /ON Pn{x) 



k=0 



r((c+l)/3m(c + 2)/3) 



Because vr„(x) satisfies (I3.29P with p = 3, b = 27, u = c + 4/3 and w = c — 1/3, we get 
our proposition from the above formula. □ 

Again, based on the above proposition, we obtain the following result in the oscillatory 
region. 

Proposition 5.4. Let x = x„ := (n + (2c + 1) /6)H with t = 108 cos^ 9,9 E [5,^-5], we 
have as n oo, 



Pn\Xn) 



3C+l/2r2(c±l)r(c±2) 



24/37r3/2(108t-t2)l/4 ^^3^^^^ 



- exp 



^TT ,^ (2c+l),^ , 2c+l, 
X cos{- - {2n + g ' )9 + (n + 



6 ^408^ ' 

(t/s)l/3 



s — 



(5.12) 



where p is defined in (I4.26p . 
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Then the Airy-type expansion can be derived. 



Theorem 5.1. Let u = n + 

respectively, we have 



2c+l 



With Kn and U{t) defined zn (1412]) and KT3\\ . 



X 



3c+ir2(£±i)r(£±2) 

4 ' e-1 

I'STT U'^ ^ t 3 



exp 



3v^r( 



-Ut-i 



6^ 

oo 



t{t - 108) 



Ai(.if/(t))5^^1i^ + Ai'(.if/(t))5: 



U 6 



Bs{U) 



(5.13) 



uniformly for < 5 < t < oo . 



Proof. When t < 108, recall the asymptotic expansion for Airy function in (12.411) and 
(I2:42|) again, we have from (jUD and fICTD 



Pn(a;) = Ci{x)Pn{x) + C2(a;)Q„(x) ~ Ci{x) 



3^6 



/2i 



3 TT 



, ^cos —i-Uit))^ - , 
[t(108-t)]i V3^ ^" 4 



- C2(a;) 



CAx) 



3v^ 



, ^sin — (-[/(t))2 - ^ 
[t(108-t)]i V3^ 4 



2v . 



3 TT 



3^6 



-C2(X 



[t(108 -t)]3 0r 
3^6 1 



cos Z/ COS 



1 A-54 



sm cos 



54 
5^ 



[t(108 -t)]4 V^r 
Comparing the above formula and (I5.12p . we have 
3^-ir2((c+l)/3)r((c + 2)/3) 



a/s(108-s) 4 



v/s(108-s) 



Ci(x) 



25/67]- 3, V 



, , a; 1 A/7rr(i) 



C2(x) = 0. (5.14) 
□ 



This completes the proof of the theorem. 

Remark 5.1. Again, we can obtain the limiting zero distribution for Q„(108ra^x) from 
the recurrence relation (11.51) . It is the same as that in (I4.29p . 

Since the above theorem is uniformly valid in the neighbourhood of the large transition 
point 108, we have the following asymptotic formula for Qn{x). 

Corollary 5.1. Let u = n + and Kn he given in (15.21) . Uniformly for a hounded real 
numher s, we have, for x = 108i^'^ + 54v^T8 s v'^ 



e-l 

X 3 exp 



^"^^^-x^ 1 QJx) 



3^r(|)" 



-iTKn- 



3--lr2(£±i)r(4^)z/i 



4^ 



TT 



Ai(s) +0(i/-3 
as u ^ 00. (5.15) 
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Proof. The proof is very similar to that in Corollary 14.11 



□ 



Remark 5.2. Again we can put the left hand side of f l5.15p into its orthonormal form. Let 
Qn{x) be the orthonormal Conrad-Flajolet polynomials II, then (I5.15P can be rewritten 

as 



e-l 

X 3 exp 



v/?r(i) . 



3^r( 



X3 QJx) 



T{c+l)u- 

23 a! 



Ai(s) + 0(z/^3) , as z/ ^ oo. (5.16) 



According to the above formula, again it seems reasonable to conjecture that at least one 
of the weight functions for the Conrad-Flajolet polynomials II should behave like 



X -i exp 



3r(i) 



-Xi 



as X — )■ oo. 



(5.17) 



Since fl5.15p is very similar to f l4.30p except some constants difference on the right 
hand side and powers of x changed on the left hand side, we have the same approximation 
about the zeros as in Corollary 14. 2[ that is, 

Xn,k = I08u^ + 54^ OfcZ/^/^ + 0(i^^/3), (5.18) 

where u = n + 

The Bessel-type expansion is also obtained as follows. 

Theorem 5.2. Let u = n + ^ . With Kn and U*{t) defined m ([ESD, flCTj) and flCTj) . 

respectively, we have 



v^r(i; 



2c+l/2p2/'c+l^-p|'c+2^ ,,S-c 



3 

Z/2 



2iV6vr 



exp 



v^r(i)_, 
3v^r( 



Z/t3 



u*{t) 



t(108-t) 



X < sm 



cos 



2lV3r(|) 3 

v^r(i) . c 

^ — ^ " ^ X 3 vr 

2iV3r(|) 3 



s=0 

oo 



s=0 

oo 



B*{U* 



s=0 



s=0 



(5.19) 



uniformly for —oo < t < M < 108. 

Proof. The proof is similar to that in Theorem 14.21 



□ 
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6 Remarks on the Moment Problem 



The general solution of an indeterminate moment problem is described in terms of four 
entire functions A{z) , B (z) , C (z) , D (z) such that 

A{z)D{z) - B{z)C{z) = 1 (6.1) 

The functions are called the Nevanlinna parametrization; see [H [19] . It is known that 
these functions have real and simple zeros. The following theorem, see [T3l Thm. 21.1.2], 
describes all solutions to an indeterminate moment problem. 

Theorem 6.1. Let JV denote the class of functions {cr}, which are analytic in the open 
upper half plane and map it into the lower half plane, and satisfy a{'z) = a{z). Then the 
formula 

df,{t-a) _A{z)-a{z)Ciz) ^ ^^2) 



z-t B{z) - a[z) D[z) 
establishes a one-to-one correspondence between the solutions fi of the moment problem 
and functions a in the class M , augmented by the constant oo. 

A theorem of Herglotz, [T9l Lemma 2.1] asserts that any a in the above theorem has 
the representation 

a{z) = —ciz + C2 + / da{t), (6.3) 

Jr z ~ i 

where a is a finite positive measure on M, ci and C2 are real constants, and ci > 0. 
The next theorem is in [3] and [T5] . 

Theorem 6.2. Let a in (16.21) be analytic inlmz > 0, and assume a maps Imz > into 
lma{z) < 0. If fi{x, a) does not have a jump at x and a{x ± iO) exist then 

dfi{x;a) a{x — iO^) — a{x -\- iO'^) 

dx ~ 2TTi \B{x) - (7{x -iO+)D{x)f' 

A consequence of Theorem 16.21 is that the polynomials associated with this moment 
are orthogonal with respect to the weight function 

'y^B^[xj + U^[x) 

for any 7 > 0. 

Berg and Pedersen |1] showed that the four entire functions A, B,C,D have the same 
order, type, and the Phragmen-Lindelof indicator. 

The formulas for the orthonormal polynomials in (I2.52p . (13.411) . (14.321) and (I5.16p . 
seem to indicate an interesting pattern about the powers of u (a constant plus n, the 
degree of the polynomial). Indeed it suggested to us to formulate the following conjecture. 
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Conjecture 6.1. Let {Pk{x)}^^-^ be a sequence of polynomials orthonormal with respect 
to the weight function w{x) on an unbounded interval andtn be the large transition point 
of Pn{x). Without loss of generality, we may further assume the interval to be (0, oo) or 
(—00,00). If the weight function w{x) has the following behavior 

— logw{x) = ©(x™) as X 00, (6.6) 

then we have 

tn = 0{n"^), as n 00. (6.7) 

2 

Moreover, uniformly for a bounded real number s, we have, for x = t„(l + sn~3) 

m;(x)5P„(x) ~ c(s) n''(l + 0(1)), as n 00, (6.8) 
where c{s) is a uniformly bounded function in s and the constant k is given by 

Here m > but is not necessarily an integer. 

Note that the above conjecture is true for Laguerre polynomials, Hermite polynomials 
and polynomials orthogonal with respect to the Freud weight e"'^'". For the asymptotics 
of Laguerre polynomials and Hermite polynomials, see [201 Sec. 8.22]. For the asymp- 
totics of polynomials orthogonal with respect to the Freud weight, see [161 Eq. (1.19)]. 

When the corresponding moment problem is indeterminate, m in (16.61) equals the 
order of any of the entire functions A,B,C or D. This seems to be the case for the 
weight function (16. 5p . It also agrees with the known weight functions for the Chen- 
Ismail polynomials and the Berg-Valent polynomials. 

In the rest of this section, we shall show that the moment problems corresponding to 
Conrad-Flajolet polynomials I and II do not have unique solutions, that is the moment 
problems are indeterminate. The following proposition is in [1] or [T9] . 

Proposition 6.1. Let P„(x) be the orthonormal polynomials satisfying the following 
recurrence relation: 

xVnix) = bn+lVn+lix) + a„P„(x) + hnVn-l{x). (6.10) 

// the moment problem has a unique solution then the series 



n=l 



diverges for z ^ M. // the above series diverges at one 2 ^ M, then the moment problem 
has a unique solution. 
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Recall that a birth and death process leads to polynomials Qn{x) defined in (11. ip . 
We have a„ = A„ + 6„ = in (16.101) and the corresponding orthonormal one 

is 



Qnix) 



-irQnix) 



n 



(6.12) 



For the orthonormal Conrad-Flajolet polynomials I, we obtain the following result from 
Theorem 14.21 

Corollary 6.1. Letu = n + ^. With and U*{t) defined m (g^D, (1419|) and (ICTj) . 

respectively, we have 



X < sm 



— cos 



v^r(|) . 

-7^ X'-i TT 

2iv^r(|) 3 



7rr(i) 1 c 

X3 TT 



3 ^'-^ V 3 

_ ;„ 2c— 1 

25/% t— 



exp 



v^r(|) 
3^r(i 



f/*(t) 



t(108-t) 



J.(z.f/*^(t))5^ 



i:(f/*) 



2.v^r(|) 



s=0 

oo 



s=0 

oo 



s=0 



(6.13) 



uniformly for — oo < t < M < 108. i/ere A„ and /i„ are defined in (II. 4p . 

Let z = z/^t be a fixed complex number, then t = zjv^ = 0{n~^). It follows from 
(I4.23P that U*{t) = 0{t) = 0{n~^). Recall the asymptotic formulas of Bessel functions 
at small arguments, we obtain 

Ji{uU*Ht)) = 0(n"i/6), J4(z/f/*^(t)) = 0(n-2/3), 
Yi{iyU*kt)) = 0(^1/^), y4(z/f/*5(t)) = 0(^2/3). 

Since A*{U*) = 1 and B*{U*) = by (I4.18p . we observe that the expressions involving 
Bessel functions and trigonometric functions are of order 0{n^^^). Next, we obtain from 
the following asymptotic estimate and (II. 4p . (14. 2 p and Striling's formula. 



Oln-"'"), K„ = 0( 



n 



-2/3N 



A combination of the above estimates yields Qn{z) = 0{n~^^^). This establishes the 
convergence of the series in (16. lip . Thus, the corresponding moment problem is indeter- 
minate. 
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For the orthonormal Conrad-Flajolet polynomials II, we obtain the following result 
from Theorem 15.21 

Corollary 6.2. Let v = n+^. With Kn and U*{t) defined in fICTD and flCTj) . 

respectively, we have 



--{-I 



X < sm 



cos 



2lv/3r(|) 3 



3 

Z/2' 



exp 



3^r(| 



-Ut3 



U*{t) 



t(108-t) 



i c 
2iV3m 3 



s=0 

oo 



s=0 



(6.14) 



uniformly for — oo < t < M < 108. i/ere A„ and /i„ are defined in (11. 5p . 



A similar argument as in the case of Conrad-Flajolet polynomials I, we conclude that 
the moment problem for Conrad-Flajolet polynomials II is also indeterminate. 
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